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ON THE DETERMINATION OF CERTAIN BASIC TYPES OF 
SUPERSONIC FLOW FIELDS*- 
By Carlo Ferrari 

SUMMARY 


A quite ^anlversal mode of attack on problems which arise In super- 
sonic flow, whether connected with flow over wings or over bodies of 
revolution, is explained, first, in great generality, and then in more 
detail, as specific applications to concrete cases are illustrated. The 
method depends on the iise of Foinrier series in the formal definition of 
the potential governing the flow and in the setting up of the boundary 
conditions. This new formulation of the many problems met in supersonic 
flow is really an exlenslon of the doublet type of "fundamental solution" 
to higher order types of singularity. The limitations and, in contrast, 
the wide field of applicability of such a means of handling these prob- 
lems with complex boundary conditions is discussed in some detail, and 
a specific example of a wing-body in-berference problem is cited as proof 
of the versatility of the method, because the results obtained by applying 
the techniques expounded herein agree well with experimentally determined 
data, even for the q^uite complex configuration t^sod to exemplify the kind 
of problem amenable to such treatment. 


1, INTRODUCTION 


For purposes of analytic trea-tment of the flow problem to be con- 
sidered here the usual rectangiilar Cartesian coordinate system is employed 
with the x-axis taken to lie in the direction of and ha-ving the same sense 

— ^ 

as the uniform (undisturbed) free-stream velocity, Voo- This free -stream 
— ^ 

velocity, Voo^ is taken to be supersonic in the discussion that follows; 
i.e., Voo > Coo where Coo denotes the velocity of sound in the undis- 
turbed stream. The flow of the gaseous fliild to be Investigated is to 

■5f* .. 

Siilla determlnazione dl alcunl tipi di carpi di corrente Iper- 
sonora, " Rendiconti dell ’ Accademla Nazionale del Llncei, Classe di 
Sclenze flslche, matematlche e natural!, serie VIII, vol. VII, no. 6; 
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be considered as resulting from the superposition upon the free-stream 
velocity of a nonuniform flow, having velocity components that are des- 
ignated as Vx, Vy, and Vg, and lying in the direction of the respec- 
tive axes (x,y,z) of the coordinate system. This nonuniform super- 
imposed flow is supposed to be small enough, in comparison with the 
speed of sound, Coo^ that it is permissible to neglect the ratios Vx/Coo^ 
Vjj/Coo, etc. in the equations governing the flow. 

It is taken for granted that, under the conditions stated above, 
there exists a velocity potential describing the flow in question, and 
in practically all cases which are of any interest for actual designs 
it will really be true that this assung)tion can be made legitimately. 

If it is then agreed that the nonuniform superimposed part of the flow 
is to be denoted by the potential. 0, it will be recognized that this 
potential will have to satisfy the relationship: 






= 0 


( 1 ) 


where the free-stream Mach niomber, I^, is defined as and, 

of course, the potential 0 will also have to obey the boundary condi- 
tions which are peculiar to each stated problem. 

A way of handling the determination of the function 0, so that it 
will satisfy equation (1) and so that it will obey the imposed boundary 
conditions, will now be explained, and its usefulness illustrated by 
consideration of problems which can be attacked by this means, both in 
the case of lifting surfaces (that is, wings) as well as in the case of 
bodies of revolution. The proposed method is based on the use of Fourier 
series. Although this technique does not afford complete universality 
in treatment of all the posed problems, as will be more clearly pointed 
out in what follows, it can be lised to fine advantage in a goodly number 
of situations by replacing the procedures which are based on the Fourier 
or the Laplace transforms (which, for that matter, have just as restricted 
limits of applicability as the analogous ones which arise in connection 
with the approach being discussed herein) or by being substituted in 
place of the techniques which stem from use of the "fundamental" (source, 
sink, doublet) solutions to equation (l), or from use of transformations 
carried out in the -complex plane. 


2. PR0BLE3® HAVING TO DO WITH FLOW OVER WINGS 

As us\ial, the wings are lma.gined to be very slender and so placed 
that the wing span lies along the y-axlsj i.e., the long dimension is 
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out the y-axis (see fig. 1 ). Let the equations which define the ventral 
and dorsal surfaces of the wing surface be, in fact, given in the form; 


z = z^(x,y) and z = Zj3^(x,y) 


and then the slenderness of the wing is supposed to be slight enough 
that the above -defined values of z will be so small at all locations 
on these surfaces as to make it possible to accept the fact that the 

derivative ^ is, to all Intents and pvirposes, equal to the direction 
ox 

cosine, with respect to the free-stream x-axis, of the normal to the 
surface. It is fiarther assumed that the wing is immersed in a stream 
of supersonic flow which has a constant value for its component lying 
in the direction of the x-axis, of magnitude Vco* The component in the 
direction of the z-axis, meanwhile, is asstmed to be known, but of rela- 
tively small size in comparison with the Voo velocity, and it may take 
on varioiis values, which will be denoted by V^'. Ifj now, the potential 

describing the flow perturbed by the wing is denoted by 0 this poten- 
tial will have to satisfy equation (l), and it will also have to conform 
to the conditions which are in^iosed at the boundaries. These further 
(boundary) conditions may be stated as follows: 

(2) Upstream of a certain sxrrface, which may be immediately defined 
just as soon as the wing-like body is specified which is to invest the 
Impinging stream, the value of 0 is zero; i.e., the basic condition is 


0 = 0 


(2) 


(3) On the wing siorface, it must be true that 

(Iz) ~ " S(x,y) ( 3 ) 

wherein the value of z to be employed is either the Zy or quanti- 

ties, depending on whether one is concerned with a point which is lying 
on the under ventral surface or on the upper dorsal svirface, respectively. 

The notation cos (n,z) signifies the cosine of the angle between the 
z-axis and the unit vector taken in the direction of the exterior normal 
to the wing surface in qviestion; i.e., this vector is represented by the 
vector n, and vinder the present hypothesis cos (n,z) = ± 1 . 
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It is convenient to dlstingiilsh between two basic types of problem 
which ccaae under this kind of analysis, and to make the differentiation 
on the basis of the sort of boundary conditions met with in each type; 
that is. 

Symmetric Types of Conflg^Jration 

In this case, the boimdary conditions to be satisfied on the wing 
may be expressed in the form 

(M') =H(x,y) (5') 

and 



-H(x,y) 


Asymmetric Types of Configiiratlon 

In this case, the boimdary conditions are expressed as 



G{x,y) 


(5") 


The first type of problem corresponds to a configuration for which 
the wing has a zero angle of attack with respect to the free-stream 

undisturbed flow, Vo^, and which possesses a symmetric profile. The sec- 
ond type of problem corresponds to a configuration for which the wing 
is a flat plate, but which has any local angle of attack whatsoever, 

with respect to the free-stream vector, Vco, so long as it is small. 


3. DEVELOPMEOW OF THE CASE OP THE SYMMETRIC TYRE OF CONFIGURATION 


In this case it will siiffice to examine the flow solely in the 

upper half -plane, where z > 0. If ,y, z) stands for the flow 

which takes place in this upper region, and if ^(2)(x,y,z) represents 
the flow in the nether region, then, of course, 


0^^^(x,y,z) = ^^^\x,y,-z) 
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The boimdaxy conditions in this case are composed of eq.iaations (5')^ 
together with the restriction that 


(for locations lying beyond the region 
occupied by the wing surface) (2') 


Now let the definition of the fimction describing the velocity com- 
ponent at the wing surface^ and also the potential fimction itself, be 
cast into the convenient forms 

H*(x,y) = Voo ^ cos = Vco ^ ^( 1 ) cos I mri 



and 


m 


wherein | = x/b, ti = y/b, and ^ = z/b, while b is a suitable length 
•used for purposes of nondlmensionadization. The value vised for b will 
be equal to the semispan of the wing in the case where the leading edge 
of the wing is supersonic everywhere, and pro'vlded that the wing tips 
are cut off in such a way that the -wing surface remains outside of the 
tip Mach cones emanating from either one of the 'wing-tlp extremities out 
at the farthest reaches of the -wing span. The value used for b will 
be larger than this semlspan Just defined, if, in contrast, these geo- 
metrical. relationships do not hold; the magnitude employed for b in 
this latter case is Illustrated in figure 2. 

Finally, it should be observed that B* is a periodic function 
of y, which is equal to the values taken on by the function H at the 
■wing's surface and it is zero for points lying out of this region, and 
this definition is to hold throughout the spanwlse interval for which 
-b < y < b. 

The fact that it is possible to -write H*'(x,y) in the foim given 
as equation (4) (i.e., the possibility of expressing the ccmponent- 
veloclty field describing the normal velocities to the ■wing surface by 
means of a Fourier series instead of in terms of a Fourier Integral) 
stems frcan the property already noted to the effect that the perturba- 
tions, which are created at any arbitrary point P(x,y) whatsoever, do 
not make themselves felt anywhere outside of the Mach cone emanating 
from P. As a result of this situation, therefore, as far as the 
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detennlnation of the field of flow about the given wing is concerned, it 
makes no difference to this flow whether one considers the wing to be 
operating by itself as an isolated entity within the impinging stream or 
whether, instead, one imagines it to be accompanied by an infinite n\jmber 
of reflections of this primary wing in the planes y = ±mb. 

If one now Inserts the second of the expressions given as equa- 
tion ( 4 ) into equation (l), it will be seen that this differential eqiia- 
tion reduces to 




- b2 


a|2 


= 


wherein = 1^2 _ 2. and where 


k replaces the constant 



( 5 ) 


Meanwhile, it is also evident that, on the basis of the first of 
the formal developments given as equation ( 4 ), the boundary condition 
reduces to 



(6) 


The expression given as eqtiation (5) above is formally analogous to 
the so-called "telegraph equation," and its solution, which is suitable 
for applying the type of boundary condition exemplified in equation (6), 
is 


m B 



bm(r) Jo 


I 


d| 


( 7 ) 


where Jq is the cylindrical Bessel function of zeroth order. 


Consequently, the vertical derivative turns out to be 


f Ji - bV 


d 5 ' 




SC 


rthjjjCl-BC) + rtk^ 


0 
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and, "becauBe of the 'boundary condition (6), it follows that 



so that the sought potential must have the form 

= %(r)jodi* (7’) 

4 . DEVELOEMENT OF THE CASE OP 1HE AS'IfMMETRIC TYPE OF CONFIGURATION 


The possibility of being able to find solutions to such asymmetric 
problems by means of the method being propounded here is restricted in 
this case to those configurations for which the leading edge as well as 
the trailing edge of the wing are supersonic, and where the wing tips 
are cut off in such a way that the wing sxarface lies outside of the Mach 
cone emanating from the very tip of the leading edge where the TnnviTmim 
span occurs . 

Under these circumstances the boundary conditions are constituted 
frcm the restrictions given as equations ( 3 ")> an.<3- of equation (2’) 
once again. If one then follows the same procedure as was utilized in 
section 3 > i't follows that the expression for the so\ight potential is 
formally given as (ref. 1) 


±ithm|7 - b|^+ 1/(1 - i'f - 


dg » 


( 8 ) 


where hj^ is, a priori, an xmdetermined function, and where it should 

be recognized that the + sign is to be employed for- the lower half -plane 
where ^ < 0, and where the - sign is to be employed for the upper half- 
plane where ^ > 0. It is evident, therefore, that the derivative of 

with respect to ^ will, be continuous eLLong the plane ^ = 0, but the 
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derivative of ^ with respect to | will be discontinuous, and the 
"Jump" will be of such size that 

/^=o" 

holds true. 

It is clearly permissible here again to concentrate attention solely 
upon the dlstiirbed flow in the upper half -plane where ^ > 0 , therefore, 
because the observation Just made above will tell one how to ccxiipute what 
the flow will be in the other lower half -plane, once the former is 
obtained. 

The boimdary conditions in this instance may now be recast into the 

form 

itBhni(0 + 1^(5') Jo ■' i - * ) - 

rtk hm(l’) Ji fk 


provided, as in the previoiis section, one sets up the convenient con- 
vention that G*'(x,y) is to represent a periodic function in y that 
is to be eqiial to the values taken on by the function c(x,y) at the 
wing's surface, and it is to be zero for points lying out of this region. 
This definition is to hold throughout the spanwise interval for which 
-b < y < b. In addition, the form of G*<-(x,y) is to be assumed, specif- 
ically, to have the appearance 


G^(x,y) = V„ G^(|) cos ^ mt) 
m 

while it has also been assumed that the derivative of a function by the 
sole parameter upon which it depends is to be denoted by a dot over the 
function, that is. 
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The integro-differeiitial eqtiation defining may also be immedi- 

ately simplified to the compressed expression 


+ It bia(l') (jo + = G&id) 


(9) 


Now apply a Laplace transformation to this integro-differential 

equation (i.e., multiply through by the factor e"^^^ and Integrate 
from 0 to 00 ^ . Thus, one obtains 


,2 ^ 


rtBhjjjP + It ^ hjjj 


/ 


. k2 


..,_i _..+s 

i'' 

W 


/ 2 ^ 

k2 


k2 

P 

i^' 

+ 

i? 


= 


where a bar over a symbol serves to indicate that this quantity stands 
for the Laplace transfoim of the function so designated. 


Standard tables of Laplace transforms could be consulted to check 
these resialts, which may now be simplified by noting that 




k2 


+ 


k2 

i? 


+ 




2p i/p*^ + 


k2 

F 



+ 


b2 



Thus the Laplace transforms of equation ( 9 ) simplifies to 


+ ^ - pj bm = Gm 

or the explicit expression for the Laplace transform of the unknown func- 
tion hja is given in the form 


itBptta + 
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so that finally one may Invert the transformation to obtain 



Once having obtained the value of hj^, it is easy to write down the 
expression for the coraponent of velocity lying in the x-dlreetion and 
located at the wing-surface, because one has simply that this component 
is given by the partial derivative of the potential 0, taken with respect 
to g, and evaluated at the plane of the wing; l.e., one finds that 

(IL ■ I"” 

T^IrthermDre, the foimiila giving the lift on the wing is just 

L = (.1)^ (11) 

m 


for m = 1,5,5,. • etc. 

where the symbol, Z, is used to denote the distance along the x-axis 
measured back from the leading edge of the root chord to the projection 
into the plane of symmetry of the trailing edge of the wing -tip profile. 

In regard to the moment taken about the y-axls, it is apparent that 
it may be computed from the relation; 

My = 2;rp^V^^b5 ^ cos ^sm dt] |hm(l)d| 


— r rl/b 


for m = 1,5,5^ • • etc. 
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5. PROBLEMS HAVUTG TO DO WITH FLOW PAST BODIES OF REVOLUTION 


The procedure discussed in the preceding sections can be extended at 
once to apply also to the solution of problems which are concerned with 
the flow over bodies of revolution. 


For this purpose let a cylindrical coordinate system (x,Y, 9 ) be 
set up, and then the equation which governs the potential, being 
sought will have the form 


5^0 

5 ? 


+ 





(15) 


Now let R = R(x) be the equation of the meridian line of the body, 
and let it be assumed that R is sxifficiently small at T locations 
along the body so that the direction cosine of the normal to this merid- 
ian line, measiired from the Y-axis, may be taken to be equal to unity. 

Furthermore, let Vnf ' stand for the ccmiponent, taken in a direc- 
tion perpendicular to the circular cross-section of the body, which 
arises from the impinging flow which invests the body. Then the boundary 
condition which must be satisfied at the surface of the body may be 
expressed mathematically by the relation 





I 


For sake of simplicity, it is also now assumed that the treatment 
to be developed is to be restricted to the case where symmetry with 

respect to the semiplanes 0 exists in the incident flow. Under 

this hypothesis it is convenient to write the normal velocity components 
and the potential being sought in the following explicit formulations: 




= 


m 


\R 


n 


sin m 0 




^ ^(x,Y)y“ Bln m 0 


( 14 ) 
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If one now inserts the second of the expressions given as equa- 
tion (l 4 ) into the differential equation governing the flow (l 3 )f li 
will he found that the defining eqimtion for the potential will have 
the form 


^ . Pm -I- 1 

8y2 Y Sr 


(15) 


and the hovmdary condition tiims out to he 



(15’) 


A suitable solution to equation ( 15 )^ which can he made to satisfy 
the houndary condition being imposed as equation (I5 ' ) , will be found 
to he 


5^m 



(16) 


where 



IqCx-BY cosh u)du 



MCA TM 1581 


13 


Thiis, the successive Individual potentials are given hy the expressions^ 


pO 

Y 01 = -B / 

Q V» 




arc cosh ^ 

l5x 


fl(x-BY cosh u) cosh u du 


pO 

Y02 “ / f2(x-BY cosh u) cosh^ u du 

"^axc cosh 

BY 


etc. 


> (IT) 


Upon in^josition of the requirement that the boundary condition ( 15 ' ) 
is to be satisfied, one obtains a set of integral equations which serve 


translator's note: It was pointed out on page 63O of an article 

by R. H. Cramer in the Journal of the Aeronautical Sciences, vol. I8, 
no. S, September 1951 > entitled "Interference Between Wing and Body at 
Supersonic Speeds - Theoretical and Experimental Determination of 
Pressures on the Body," that the res^ilt given here for for m > 1 , 

is incorrect j the correct formula is, for m = 2, 


Y^02 = / f2(x-BY cosh u) ^2 cosh^ u - l) du 

^arc cosh ^ 

while, in general, the use of hyperbolic functions of multiples of the 
argtmient u gives a more ccmpact form, which is easy to work with; i.e., 
in general it is true that 



fj^(x-BY cosh u) (cosh m u) 


du 
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to determine the arhitrary functions which are a priori unknown. 

Thus, applying these conditions, one finds that^ 


B- 


iHH-l 


(_i)im-i 

O V, 


arc cosh ^ 


f_i(x-BR cosh u) 


(cosh du = -Pjjj (18) 


The determination of the values of the fj^'s appearing in formula (18) 

may be carried out by using a step-by-step procedure which is entirely 
analogous to the one employed by Von KdrmAn in his work on determining 
the flow about a body of revolution at zero angle of attack. 

It is Important to point out that if one only has in mind to calcu- 
late the force distribution along the axis of the body and the corre- 
sponding moment, and if one is not interested in knowing the local veloci- 
ties or pressiires aroimd the body, then it is merely necessary to 
calculate </)q and 02* 


6. APPLICATIONS 


The procedtore that has been prc^omded above has been applied (ref. 2) 
to the situation arising in the study of the question of wing-body inter- 
ference. The wing -body configiiratlon considered in this particular appli- 
cation of the method is depicted in the appended figure 5* 

The wing used in this configuration is a flat plate, whose semispan 
is equal to 4 Rq, where Rq Is the raditis of the circular cross-section 

taken through the body at the location where the body is widest. The 
leading edge of this wing is located 5^0 downstream from the tip of the 
nose of the body. The free-stream flow is inrplnglng on the body at a 
speed which is twice the speed of somd in the undisturbed stream. 

^Translator 's note: In view of the correction pointed out in Note 1 

above, it will be seen that this formula for determining the functions 

is also Incorrect, except for m = 1; for higher Integral values of m, 
the correct formula is : 


B“^^ (-1)^4 / fjj^(x-BR cosh u) [cosh m u cosh du = -Fj^ 

'-'arc cosh 
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Tlie ciirves stLOwn in figure 5 give the value of the pressure coeffi- 
cient, ~ out along the span of the wing, in the mid-chord loca- 

tion (i.e., along the wing axis), for points on the upper (dorsal) side 
of the wing. These coefficients have been calculated by the method out- 
lined in section h, and there are shown results for various angles of 
attack, which apply to such points on the upper side of the wing pro- 
files at their mid-chord positions. 

In addition, some ejjperimental test points obtained by R. H. Cramer 
(see ref. 2 ) are also plotted on these curves. These results were 
obtained from experiments carried out in the supersonic tunnel of the 
Dalngerfleld Aeronautical Laboratory. 

The agreement between the computed and experimentally determined 
resiilts is very good from a qualitative viewpoint. In tegard to the 
more precise details of the quantitative comparison between the resTolts 
it is worthy of note that the experimental results exhibit a certain 
amount of dissymmetry as one passes frcm positive angles of attack to 
negative angles of attack. Such a dissymmetry cannot be predicted, or 
should not be expected, from the type of theoretical treatment being 
considered here. 

In order to bring about a more valid comparison of these results, 
it would appear logical, in face of such evident dissymmetry, to take 
for the representative eipierlmental value, at a given value of the angle 
of attack, 3, the one which is obtained from averaging the result 
obtained for an angle of attack equal to +3 with the result obtained 
at -3. Such average values have been computed and are designated in 
the plots of figure 5 means of solid circles. These adjusted values 
lie much closer to the theoretically derived curves at almost all 
locations . 


Translated by R. H. Cramer 

Cornell Aeronautical Laboratory, 3 nc. 

Buffalo, New York 
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Figure 1.- Orientation of coordinate axes and location of typical wing 

plan form therein. 



5’igure 2. - Definition of the interval of periodicity required for application 
of the Fourier series technique when leading edges are subsonic. 
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Figure 3 .- Pressure distribution along the wing axis: Comparison of 
experimental results with predictions based on the method expounded 
in section 4 . 
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